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We report on a search for the production of new resonance states in the process 77 — > DD. A 
candidate C-even charmonium state is observed in the vicinity of 3.93 GeV /(? . The production rate 
and the angular distribution in the 77 ccntor-of-mass frame suggest that this state is the previously 

unobserved Xc2) the 2^ P2 charmonium state. 

PACS numbers; 13.66.Bc, 14.40.Gx 



INTRODUCTION 

The masses and other properties of the radial-ground 

and radially excited states of charmonium provide valu- 
able input to QCD models that describe heavy quarko- 



nium systems. To date, radial excitation states of char- 
monium have been found only for the "^^^^Lj = ^Si 
states {tp) and ^Sq states {rjc)- No radially excited ^Pj 
states (xcj) have yet been found, even though the three 
radial ground states have been already well established. 
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The first radially excited Xcj states are predicted to 
have masses between 3.9 and 4.0 GeV/c^ [1, 2], which 
is considerably above DD threshold. If the masses of 
these states lie between the DD and D*D* thresholds, 
the Xco(2P) ix'co) and_Xc2(2F) (Xc2) are expected to de- 
cay primarily into DD, although the x'c2 could also de- 
cay to DD* if it is energetically allowed. (The inclusion 
of charge-conjugated reactions is implied throughout this 
paper.) Recently, two new charmonium-like states in this 
mass region, the X(3940) [3] and y(3940) [4], were re- 
ported by Belle. Decays of either of these states to DD 
have not been observed [3] . 

In this paper we report on a search for the x'cj { J = 
or 2) states and other C-even charmonium states in the 
mass range of 3.73 - 4.3 GeV/c^ produced via the 77 — > 
DD process. The results presented here are preliminary. 



Kaon candidates arc separated from pions based on 
normalized kaon and pion likelihood functions obtained 
from the particle identification system (Lk and L,r, re- 
spectively) with a criterion, Lk/(Lk + ^tt) > 0.8. All 
tracks that are not identified as kaons are treated as pi- 
ons. 

Signal candidates are triggered by a variety of track- 
triggers that require two or more CDC tracks with as- 
sociated TOF hits, ECL clusters or a minimum sum of 
energy in the ECL. For the four and six charged track 
topologies used in this analysis, the trigger conditions 
are complementary to each other and, in combination, 
provide a high trigger efficiency 95%). 



EVENT SELECTION 



DATA AND DETECTOR 

The analysis uses data recorded in the Belle detector at 
the KEKB asymmetric collider [5]. The data sam- 

ple corresponds to an integrated luminosity of 280 fb~^, 
accumulated on the T(45) resonance {^/s = 10.58 GeV) 
and 60 McV below the resonance. Since the beam en- 
ergy dependence of two-photon processes is small, we 
combine both samples. We study the two-photon pro- 
cess e'^e~ e~^e~DD in the "no-tag" mode, i.e. where 
neither the final-state electron nor positron is detected. 
We restrict the virtuality of the incident photons to be 
small by imposing a strict requirement on the transverse- 
momentum balance of the final-state hadronic system 
with respect to the beam axis . 

A comprehensive description of the Belle detector is 
given elsewhere [6]. We mention here only the detector 
components essential for the present measurement. 

Charged tracks arc reconstructed from hit information 
in a central drift chamber (CDC) located in a uniform 
1.5 T solenoidal magnetic field. The z axis of the detector 
and the solenoid are along the positron beam, with the 
positrons moving in the —z direction. The CDC mea- 
sures the longitudinal and transverse momentum com- 
ponents (along the z axis and in the np plane, respec- 
tively). Track trajectory coordinates near the collision 
point are measured by a silicon vertex detector (SVD). 
Photon detection and energy measurements are provided 
by a CsI(Tl) electromagnetic calorimeter (ECL) . Species 
of charged hadron are identified by means of informa- 
tion from timc-of-flight counters (TOF) and a silica- 
aerogel Cherenkov counters (ACC). The ACC provides 
separation between kaons and pions for momenta above 
1.2 GeV/c. The TOF system consists of a barrel of 128 
plastic scintillation counters, and is effective for K/tt 
separation for tracks with momenta below 1.2 GeV/c. 
Low energy kaons are also identified by specific ioniza- 
tion (dE/dx) measurements in the CDC. 



We search for exclusive DD production in the following 
four combinations of decays: 

77 -» D°D°, D° K-TT+, D° K+n- (N4), 
77 ^ D°D°, D° K-TT+, £>° ^ K+TT-TT^ (N5), 
77 ^ D°D°, D° K-n+, £>° ^ K+tt-tt+tt- (N6), 
77 ^ D+D-, D+ K-n+n+, D' K+n'-K- (C6). 

The symbols in parentheses are used to designate each 

of the final states. For the four-prong processes (N4 and 
N5) the selection criteria are: four charged tracks, each 
one with (L) a transverse momentum in the laboratory 
frame of pt > 0.1 GcV/c and a distance of closest ap- 
proach to the nominal collision point of |dr| < 5 cm and 
\dz\ < 5 cm; the absolute value of the average dz for 
all of the tracks, \dz\ < 3 cm; two or more of the four 
tracks must have (S) pt > 0.4 GeV/c, \dr\ < 1 cm, and 
-0.8660 < cos6i < -F0.9563, where 6 is the laboratory 
frame polar angle; no photon clusters with an energy 
greater than 400 MeV; the charged track system consists 
of a K~^K~n~^'jT~ combination; the larger of the two neu- 
tral Kn invariant masses should lie within ±15 MeV/c^ 
of the nominal mass. For the N4 process, we require 
that the smaller neutral Ktt mass is within til MeV/c^ 
of the nominal mass. For the N5 process, we re- 
quire that the remaining Ktt combination, when com- 
bined with a TT^ candidate, has an invariant mass in the 
range 1.83 < M(K+7r-7r") < 1.89 GcV/c^. A tt° can- 
didate is any pair of photons in the event that fits the 
— > -y-y hypothesis with < 4. If there are multiple tt" 
candidates, we select the one that results in M(iir+7r~7r°) 
closest to the nominal mass. 

For the six- prong processes (N6 and C6), we re- 
quire exactly six tracks with particle assignments 
i4r"'"if~7r+7r~7r+7r~, where all six pass the looser and two 
to four pass the more stringent track criteria that are in- 
dicated by (L) and (S) above, respectively. For the N6 
process, one combination is required to have |AM|i = 
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|M(K'+7r") - m£,o| < 15 MeV/c^ while the remain- 
ing tracks have |AM|2 = \M {K~ tt^ i:~ -n^) - m^ioj < 
30 MeV/c^. When there are multiple combinations, we 
take the one with the smallest |AM|i + |AAf|2. For 
the C6 process, we require \M{K^-k'^tt^) — m^j+j < 
30 MeV/c^ for each of the charge combinations, where 
m£)+ is the nominal -D"*" mass. 

The invariant mass distributions for D-meson candi- 
dates reconstructed according to the above selection pro- 
cedure are shown in Fig. 1. For all processes, we require 
that there are no extra 7r° candidates with transverse 
momenta larger than 100 MeV/c. We also apply the 
following kinematical requirements to reject initial-state 
radiation or pseudo-Compton processes (ISR veto). The 
invariant mass constructed from all of the tracks accepted 
by the more stringent criteria is less than 4.5 GeV/c^ 
(here a zero rest mass is assigned to each charged track), 
and the missing mass squared of the system recoiling 
against the detected tracks is larger than 2 (GeV/c^)^. 
The DD candidate system is also required to satisfy: 
P,{DD) > {M{DDf - 49 GeVVc'^)/(14 GcV/c^) + 
0.6 GeV/c, where Pz {DD) and M{DD) arc the momen- 
tum component in the z direction in the laboratory frame 
and the invariant mass, respectively. This condition elim- 
inates the ISR events from e+e" D^*^D'^*^^ efficiently, 
in case the photon is emitted in the forward direction 
with respect to the incident electron. We compute the 
invariant mass of the DD system using the measured 
3- momenta of each D candidate (Pd) and energy deter- 
mined from Ed = \J + "t-ij , where niD is the nominal 
mass of the neutral or charged D. 

We calculate Pt{DD), the total transverse momentum 
in the e^e~ center-of-mass (cm.) frame with respect to 
the incident e^e" axis that approximates the direction 
of the two-photon collision axis. In the two-dimensional 
region M{DD) < 4.3 GeV/c^ and Pt{DD) < 0.2 GeV/c, 
we find 159 N4-process events, 110 N5-process events, 
240 N6-process events and 86 C6-process events. 

DISTRIBUTIONS OF THE DD CANDIDATES 

In Fig. 2 we show the M{DD) distributions separately 
for D°L>o (sum of N4, N5 and N6) (Fig. 2(a)) and D+D~ 
(Fig. 2(b)) and for the combined charged and neutral 
channels (Fig. 2(c)). In the figures, two event concen- 
trations are evident: one near 3.80 GeV/c^ and another 
near 3.93 GeV/c^. Here, we have appHed the requirement 
Pt{DD) < 0.05 GeV/c to enhance exclusive two-photon 
77 DD production. 

The invariant-mass distribution for the combined 
D°D° and D^D~ channels is shown for 10-MeV/c^ bins 
in Fig. 3. The curve is the result of an unbinned like- 
lihood fit to the data in the region 3.80 < M{DD) < 
4.10 GeV/c^ using a relativistic Breit-Wigner signal 
function plus a background of the form ^ M{DD)^°' , 
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FIG. 1: (a) Invariant mass distribution of K^tt^ for N4 can- 
didate events, (b) Invariant mass distribution of K^tv^tt''' in 
N5 candidate events where an accompanying TT 
(tagged D) is found. (c) Invariant mass distribution of 
K^Tv^ in N6 candidates, (d) Invariant mass distribution of 
K'^ -K~ ■k'^ -K~ in N6 candidates with a tagged D" K~tx^ . 
(e) Invariant mass distribution of _ft'^7r~7r~ in C6 candidates 
with a tagged K^tt^tt'^ . 

where a is a free parameter. The mass dependence 
of the efficiency and the two-photon luminosity func- 
tion is taken into account in the fit. These are com- 
puted using the TREPS Monte-Carlo(MC) program [7] 
for e+e~ e'^e~DD production together with JET- 
SETT. 3 decay routines [8] for the D meson decays (using 
PDG2004 [9] values for the decay branching fractions). 
The M{DD) dependence of this product is not large. (At 
M (DD) = 3.93 GeV/c^, there is a 13% decrease in this 
product for a 0.1 GeV/c^ increase in M{DD).) 

The results of the fit for the resonance mass, width 
and total yield of the resonance are M = 3931 ± 
A{stat) MeV/c2, F = 20 ± d,{stat) MeV and 41 ± 
\\{stat) events, respectively. The mass resolution, which 
is estimated by MC to be 2-3 MeV/c^, is neglected in 
the fit. The statistical significance of the peak is 5.5cr, 
which is derived from the square root of the difference 
of the logarithmic-likelihoods for fits with and without a 
resonance peak component, shown in the figure as solid 
and dashed curves, respectively. 

Systematic errors for the parameters M and F are 
2 MeV/c^ and 3 MeV, respectively. The former is domi- 
nantly due to the uncertainty in the mass of the D mesons 
(1 MeV/c^ for the resonance mass) and the choice of the 
Breit-Wigner function formula (1 MeV/c^). We consider 
here several different Breit-Wigner functional forms for 
spin and 2 resonances, phase-space and wave-function 
variations. In the latter, we also consider the effects of 
the finite invariant-mass resolution in the fit. 

The Pt{DD) distribution in the peak region, 3.91 < 
M{DD) < 3.95 GeV/c^, is shown in Fig. 4. Here the 
Pt requirement has been relaxed. The experimental data 
are fitted by a shape that is expected for exclusive two- 
photon DD production plus a linear background. We 
expect non-charm and non-exclusive backgrounds to be 
nearly linear in Pt{DD). The fit uses a binned-niaximum 
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FIG. 2: Invariant mass distributions for the (a) D°D° chan- 
nels and (b) the D+D' mode, (c) The combined M{DD) 
distribution. 
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FIG. 4: The experimental Pt{DD) distribution (points with 
error bars) for events in the 3.91 < M{DD) < 3.95 GeV/c^ re- 
gion and the fit (histogram) based on the exclusive 77 — > DD 
process MC plus a linear background (dotted line). The dot- 
dashed line shows the location of the Pt selection requirement. 
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FIG. 3: The sum of the M{DD) invariant mass distributions 
for all four processes. The curves show the fits with (solid) and 
without (dashed) a resonance component. The histograms 
show the distribution of the events from the D-mass sidebands 
(see the text). 



likelihood method with the scales of the two components 
treated as free parameters. The linear-background com- 
ponent, 0.6 ±0.4 events for Pt{DD) < 0.05 GcV/c^, and 
the goodness of fit, ~x^/d.o.f = 18.7/18, indicate that 
the events in the peak region originate primarily from 
exclusive two-photon events. 

The Pt{DD) distribution produced by DD* and D*D* 
events is expected to be distorted by the transverse mo- 
mentum of the undetected slow pion(s), which peaks 
around 0.05 GeV/c. Such a distortion is not evident 
in the observed distribution. Prom a fit that includes 
a D'^*^D* component, together with DD production and 
a linear background, we find that at the 90% confidence 
level, less than 6 of the 46 events observed in the selected 
Pt region originate from D^*^D*. 

We investigate possible backgrounds from non-DD 



sources using D-sideband events. The histograms in 
Fig. 3 show the invariant mass distributions for events 
where the D-meson is replaced by a hadron system from 
a D-signal mass sideband region below the signal region 
with the same width as the signal mass region. Wc use 
two types of sideband events: one where one Z?-meson 
candidate is in the signal mass region (shaded histogram), 
and another whore both entries arc from the sidebands 
(open- histogram). Since there is no significant difference 
between the two sideband distributions, we conclude that 
the sideband events are dominated by non-charm back- 
grounds. Since the distributions of the two kinds of the 
sideband events are essentially equivalent, we combine 
them and scale by a half in order to compare to the DD 
signal yield. 

The sideband-event distribution has no enhancement 
in the peak region but does include a broad cluster 
around 3.80 GeV/c^ with a shape that is similar to the 
low mass enhancement seen for the DD candidates. Since 
the level of the sideband events is only 10-20% of the DD 
candidates, we conclude that the lower mass enhance- 
ment is dominantly DD (inclusive or exclusive) events. 

Figures 5(a) and (b) show the M{DD) distributions for 
events with | cos^^*| < 0.5 and | cos^?*| > 0.5, respectively, 
where 9* is the angle oi a, D meson relative to the beam 
axis in the 77 cm. frame. It is apparent that the events 
in the 3.93 GeV/c^ peak tend to concentrate at small 
I cos 0* I values. 

The points with error bars in Fig. 5(c) show the event 
yields in the 3.91 GeV/c^ to 3.95 GeV/c^ region ver- 
sus I cos 9* I . Background, estimated from events in the 
M{DD) sideband, is indicated by the histogram. A MC 
study indicates that the efficiency is uniform in | cos0*|. 
For a spin-0 resonance this distribution should be flat. In 
contrast, a spin-2 resonance is expected to be produced 
with helicity-2 along the incident axis [10, 11], in which 
case the expected angular distribution is oc sin^ 9* . 
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FIG. 5: M{DD) distributions for (a) Icos^*] < 0.5 and 
(b) I cos I > 0.5. (c) The |cos^*| distributions in the 
3.91 < M{DD) < 3.95 GeV/c^ region (points with error bars) 
and background scaled from the M{DD) sideband (sohd his- 
togram). The sohd and dashed curves are expected distri- 
butions for the spin=2 (helicity=2) and spin=0 hypotheses, 
respectively, plus a curve that interpolates the non-peak back- 
ground (dotted), with total area normalized to the observed 
number of events. 



The solid curve in Fig. 5(c) shows the expectation 
using sin** 9* to represent the signal plus a term pro- 
portional to 1 + a cos^ 6* that interpolates the back- 
ground (dotted curve). The comparison to the data has 

/d.o.f. ~ 7.3/9. Here the functions arc normalized to 
the total numbers of signal and background events. (We 
fix the number of the background events to be 13, which 
is obtained from the fit to the invariant-mass distribu- 
tion.) A comparison using a constant term to repre- 
sent the signal (dashed curve) gives a much poorer fit: 
X^/rf.o./. = 32.2/9. The data significantly prefer a spin 
two assignment over spin zero. 



TWO-PHOTON DECAY WIDTH 

No charmonium state that decays into DD with a mass 
near 3.93 GeV/c^ has been previously reported. We find 
no corresponding event concentration in this mass region 
in the sample of ISR events (normally rejected by the ISR 
veto), which would be expected in case of production of 
a J^'-^ = \~~ meson (■0). 

Using the number of observed signal events and the 
branching fractions and efficiencies for the four decay 
channels, we determine the product of the two-photon 
decay width and DD branching fraction (multiplied by 
the spin factor) to be (2 J-h l)r^-y(.^(3930))B(Z(3930) ^ 
DD) = 1.13 ± q.30(stai.) keV. Here, we define 
S(Z(3930) ^ DD) = S(Z(3930) ^ L>°r>°) + 
;B(Z(3930) ^ D+D-) and assume H(Z(3930) ^ 
D"D") = S(Z(3930) D+D-) according to isospin 
invariance, where Z(3930) is used as a tentative desig- 
nation for the observed state. 

The observed signals for the D"_D" and D^D~ modes 
are consistent with isospin invariance. The results on 



mass, decay angular distributions and F-^-^ [12] are all 
consistent with expectations for the x'c2^ the 2^P2 char- 
monium state. 

We assign a 16% total systematic error to the present 
measurement. This is primarily due to uncertainties 
in the track reconstruction efficiency (7%), luminosity 
function (5%), MC statistics (7%) and the D-meson 
branching fractions (9%), added in quadrature. Using 
J = 2, the above result gives V^^B{Z{'3,Q'iQ) DD) = 
0.23 ± OM{stat) ± Om{sys) keV. 



CONCLUSION 

We have observed an enhancement in DD invariant 
mass near 3.93 GeV/c^ in 77 DD events. The statis- 
tical significance of the signal is 5.5a. The observed an- 
gular distribution is consistent with two-photon produc- 
tion of a tensor meson. Preliminary results for the mass, 
width, and the product of the two-photon decay width 
times the branching fraction to DD are: M = 3931 ± 
4{stat) ± 2{sys) MeV/c^, F = 20 ±8{stat) ± 3(sys) MeV 
and T^^B{-^ DD) = 0.23 ± 0.06{stat) ± 0.04(sys) keV 
(assuming J = 2), respectively. The measured proper- 
tics arc consistent with expectations for the previously 
unseen Xc2 charmonium state. 
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